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1 Introduction 



Galaxies are attractive objects to study. The magnetism of their natural 
t^. ! beauty adds to the fascinating diversity of physical processes that occur over 

an enormous range of scales from the global dimension of order 10 kpc 1 down 
to the viscous turbulent scales of 1000 km and less. The visual image of a 
galaxy (see Fig. |T} is dominated by the optical light mostly produced by stars 
^ ■ that contribute most of the visible galactic mass (2 x 10 11 M Q for the Milky 

Way, where M & = 2 x 10 30 kg is the mass of the Sun). A few percent of the 
galactic mass is due to the interstellar gas that resides in the gravitational 
field produced by stars and dark matter. Spiral galaxies are flat (Fig. [TJ) 
because the stars and gas rapidly rotate. The gas is ionized by the UV and 
X-ray radiation and by cosmic rays; the degree of ionization of diffuse gas 
^ . ranges from 30% to 100% in various phases — see Sect. I2.1|) . Interstellar 

gas is involved in turbulent motions that can be detected because the asso- 
ciated Doppler shifts broaden spectral lines emitted by the gas beyond their 
width expected from thermal motions alone. The effective mean free path of 
interstellar gas particles is small enough to justify a fluid description under 
a broad range of conditions. Altogether, interstellar gas can be reasonably 
described as an electrically conducting, rotating, stratified turbulent fluid 
— and thus a site of MHD processes discussed elsewhere in this volume, 
including various types of dynamo action. 



1 A length unit appropriate to galaxies is lkpc w 3.1 x 10 19 m 3262 light years. The 
distance of the Sun from the centre of the Milky Way is % ~ 8.5 kpc. 
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Figure 1: Optical images of two nearby spiral galaxies. Left: M51, the 
Whirlpool galaxy (with a satellite galaxy at the top). Right: NGC 891 (Cour- 
tesy of the Canada-France-Hawaii Telescope/ J. -C. Cuillandre/Coelum). 
M51 is one of nearby galaxies (distance 9.6 Mpc) notable for its prominent 
spiral pattern. M51 is the first galaxy where a well ordered, large-scale mag- 
netic field was detected (Segalovitz et al, 1976) and studied in fine detail. 
NGC 891 is at about the same distance as M51, but seen nearly edge-on, so 
the thinness of the galactic disc is evident. The dark strip along the galactic 
disc and filaments extended away from the galactic plane are due to obscura- 
tion by interstellar dust. The filaments trace gas outflow from the disc into 
the halo. 
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The energy density of interstellar magnetic fields is observed to be com- 
parable to the kinetic energy density of interstellar turbulence and cosmic 
ray energy density, and apparently exceeds the thermal energy density of 
interstellar gas (Cox, 1990). Therefore, interstellar gas, magnetic field and 
cosmic rays form a complex, nonlinear physical systems whose behaviour is 
equally affected by each of the three components. The system is so complex 
that magnetic fields and cosmic rays — the components that are more dif- 
ficult to observe and model — are often neglected. Such a simplification is 
perhaps justifiable at very large scales of order 10 kpc, where the motions of 
interstellar gas (mainly the overall rotation) are governed by gravity: sys- 
tematic motions at a speed in excess of 10-30 km s _1 are too strong to be 
affected by interstellar magnetic fields. However, motions at smaller scales 
(comparable to and less than the turbulent scale, i ~ 0.1 kpc) are strongly 
influenced by magnetic fields. In particular, interstellar turbulence is in fact 
an MHD turbulence. In this respect, the interstellar environment does not 
differ much from stellar and planetary interiors. 

Until recently, interstellar magnetic fields had been a rather isolated area 
of galactic astrophysics. The reason for that was twofold. Firstly, magnetic 
fields are difficult to observe and model. Secondly, they were understood 
too poorly to provide useful insight into the physics of interstellar gas and 
galaxies in general. The widespread attitude of galactic astrophysicists to 
interstellar magnetic fields was succinctly described by Woltjer (1967): 

The argument in the past has frequently been a process of elim- 
ination: one observed certain phenomena, and one investigated 
what part of the phenomena could be explained; then the unex- 
plained part was taken to show the effects of the magnetic field. It 
is clear in this case that, the larger one's ignorance, the stronger 
the magnetic field. 

The attitude hardly changed in 20 subsequent years, when Cox (1990) ob- 
served that 

As usual in astrophysics, the way out of a difficulty is to invoke 
the poorly understood magnetic field. . . . One tends to ignore the 
field so long as one can get away with it. 

The situation has changed dramatically over the last 10-15 years. Theory and 
observations of galactic magnetic fields are now advanced enough to provide 
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useful constraints on the kinematics and dynamics of interstellar gas, and 
the importance and role of galactic magnetic fields are better appreciated. 

In this chapter, we review in Sect. [2]those aspects of galactic astrophysics 
that are relevant to magnetic fields, and briefly summarize in Sect. H3 our 
observational knowledge of magnetic fields in spiral galaxies. Section 0] is 
an exposition of the current ideas on the origin of galactic magnetic fields, 
including the dynamo theory. The confrontation of theory with observations 
is the subject of Sect. El where we summarize the advantages and difficulties 
of various theories and argue that the mean-field dynamo theory remains the 
best contender. Magnetic fields in elliptical galaxies are briefly discussed in 
Sect. E 

2 Interstellar medium in spiral galaxies 

2.1 Turbulence and multi-phase structure 

The interstellar medium (ISM) is much more inhomogeneous and active than 
stellar and planetary interiors. The reason for that is ongoing star formation 
where massive young stars evolve rapidly (in about 10 6 yr) and then explode 
as supernova stars (SN) releasing large amounts of energy (Esn ~ 10 51 erg 
per event). These explosions control the structure of the ISM. 

SN remnants are filled with hot, overpressured gas and first expand super- 
sonically; at this stage the gas surrounding the blast wave is not perturbed. 
However, a pressure disturbance starts propagating faster than the SN shell 
as soon as the expansion velocity becomes comparable to or lower than the 
speed of sound in the surrounding gas — at this stage the expanding SN 
remnant drives motions in the surrounding gas, and its energy is partially 
converted into the kinetic energy of the ISM. When pressure inside an SN 
remnant reduces to values comparable to that in the surrounding gas, the 
remnant disintegrates and merges with the ISM. Since SN occur at (almost) 
random times and positions, the result is a random force that drives random 
motions in the ISM that eventually become turbulent. The size of an SN 
remnant when it has reached pressure balance determines the energy-range 
turbulent scale, 

0.05-0.1 kpc. 

A useful review of supernova dynamics can be found, e.g., in Lozinskaya 
(1992), and the spectral properties of interstellar turbulence are discussed by 
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Armstrong et al. (1995). Among numerous reviews of the multi-phase ISM 
we mention that of Cox (1990) and a recent text of Dopita & Sutherland 
(2003). 

About / = 0.07 of the SN energy is converted into the ISM's kinetic 
energy. With the SN frequency of i^sn ~ (30 yr) -1 in the Milky Way (i.e., 
one SN per 30 yr), the kinetic energy supply rate per unit mass is esN = 
/^SM-^SN^g^s ~ 10" 2 ergg _1 s _1 , where M gas = 4 x 10 9 M Q is the total mass 
of gas in the galaxy. This energy supply can drive turbulent motions at a 
speed Mo such that 2wq/£ = esN (where the factor 2 allows for equal contri- 
butions of kinetic and magnetic turbulent energies), which yields 

Mo ~ 10-30 km s _1 , 

a value similar to the speed of sound at a temperature T = 10 4 K or higher. 
The corresponding turbulent diffusivity follows as 

77t ~ \£u « (0.5-3) x 10 26 cm 2 s" 1 . (1) 

Supernovae are the main source of turbulence in the ISM. Stellar winds is 
another significant source, contributing about 25% of the total energy supply 
(e.g., §VI.3 in Ruzmaikin et al., 1988). 

The time interval between supernova shocks passing through a given point 
is about (McKee & Ostriker, 1977; Cox, 1990) 

r = (0.5-5) x 10 6 yr. 

After this period of time, the velocity field at a given position completely 
renovates to become independent of its previous form. Therefore, this time 
can be identified with the correlation time of interstellar turbulence. The 
renovation time is 2-20 times shorter than the 'eddy turnover' time £/uo ~ 
10 7 yr. This means that the short-correlated (or ^-correlated) approximation, 
so important in turbulence and dynamo theory (e.g., Zeldovich et al., 1990; 
Brandenburg & Subramanian, 2004), can be quite accurate in application to 
the ISM — this is a unique feature of interstellar turbulence. Note that the 
standard estimate (JTJ) is valid if the correlation time is £/uq. If the renovation 
time was used instead, the result would be t/t ~ C 2 jr ~ 10 27 cm 2 s _1 , a value 
an order of magnitude larger than the standard estimate. 

Another important result of supernova activity is a large amount of gas 
heated to t — 10 6 K (Fig. |2J). The gas is so tenuous that the collision rate 
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Figure 2: SN remnants are expanding bubbles of hot gas that emits thermal 
X-rays. This is illustrated by the X-ray image of Tycho's supernova remnant 
(left panel; courtesy of the ROSAT Mission and the Max-Planck-Institut fur 
extraterrestrische Physik) whose parent star's explosion in 1572 was recorded 
by the famous Danish astronomer Tycho Brahe. The hot gas cools only 
slowly, and SN remnants often merge. The right panel shows a false-colour 
optical (Ha) image of two SN remnants DEM L316 in the Large Magellanic 
Cloud which appear to be colliding (Williams et al, 1997; image produced 
by the Magellanic Cloud Emission-Line Survey, reprinted with permission). 

Table 1: The multi-phase ISM. The origin and parameters of the most impor- 
tant phases of interstellar gas: n, the mid-plane number density in hydrogen 
atoms per cm 3 ; T, the temperature in K; c s , the speed of sound in kms -1 ; 
h, the scale height in kpc; and fy, the volume filling factor in the disc of the 
Milky Way, in per cent. 



Phase 


Origin 


n 


T 


c s 


h 


fv 


Warm 




0.1 


10 4 


10 


0.5 


60-80 


Hot 


Supernovae 


lO- 3 


10 6 


100 


3 


20-40 


Hydrogen clouds 


Compression 


20 


10 2 


1 


0.1 


2 


Molecular clouds 


Self-gravity, 


10 3 


10 


0.3 


0.075 


0.1 



thermal instability 
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of the gas particles is low, and so its radiative cooling time is very long 
and exceeds r: the hot bubbles produced by supernovae can merge before 
they cool (Fig. |2J. A result is a network of hot tunnels that form the hot 
component of the ISM. Altogether, the interstellar gas is found in several 
distinct states, known as 'phases' (this usage may be misleading as most of 
them are not proper thermodynamic phases) whose parameters are presented 
in Tabled Some of the parameters (especially the volume filling factors) are 
not known confidently, so estimates of Table ^ should be approached with 
healthy skepticism. The warm diffuse gas can be considered as a background 
against which the ISM dynamics evolves; this is the primary phase that 
occupies a connected (percolating) region in the disc, whereas the hot gas 
may or may not fill a connected region. The warm gas is ionized by the 
stellar ultraviolet radiation and cosmic rays; its degree of ionization is about 
30% at the Galactic midplane. The hot gas is so hot that it is fully ionized 
by gas particle collisions. 

The locations of SN stars are not entirely random: 70% of them cluster in 
regions of intense star formation (known as OB associations as they contain 
large numbers of young, bright stars of spectral classes O and B) where 
gas density is larger than on average in the galaxy. Collective energy input 
from a few tens (typically, 50) SN within a region about 0.5-1 kpc in size 
produces a superbubble that can break through the galactic disc (Tenorio- 
Tagle & Bodenheimer, 1988). This removes the hot gas into the galactic 
halo and significantly reduces its filling factor in the disc (from about 70% 
to 10-20%). This also gives rise to a systematic outflow of the hot gas to 
large heights where the gas eventually cools, condenses and returns to the 
disc after about 10 9 yr in the form of cold, dense clouds of neutral hydrogen 
(Wakker & van Woerden, 1997). This convection- type flow is known as the 
galactic fountain (Shapiro & Field, 1976), and it can plausibly support a 
mean-field dynamo of its own (Sokoloff & Shukurov, 1990). Another aspect 
of its role in galactic dynamos is discussed in Sect. 14.31 The vertical velocity 
of the hot gas at the base of the fountain flow is 100-200 km s _1 (e.g., Kahn 
& Brett, 1993; Korpi et al, 1999a,b). 

2.2 Galactic rotation 

Spiral galaxies have conspicuous flat components because they rotate rapidly 
enough. The Sun moves in the Milky Way at a velocity of about u & = 
SqQq = 220 km s -1 , to complete one orbit of a radius s ~ 8.5 kpc in 
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27r/f2 = 2.4 x 10 8 yr. These values are representative for spiral galaxies 
in general. The Rossby number is estimated as 

Ro Uo > 



The vertical distribution of the gas is controlled, to the first approximation, 
by hydrostatic equilibrium in the gravity field produced by stars and dark 
matter, with pressure comprising thermal, turbulent, magnetic and cosmic 
ray components in roughly equal proportion (e.g., Boulares & Cox, 1990; 
Fletcher & Shukurov, 2001). The semi-thickness of the warm gas layer is 
about h = 0.5 kpc, i.e., the aspect ratio of the gas disc is 

e = — ~ 0.06 . (2) 

S & 

Since the gravity force decreases with radius s together with the stellar mass 
density, h grows with s at s > 10 kpc (see §VI.2 in Ruzmaikin et al, 1988, 
for a review). 

However, the hot gas has larger speed of sound and turbulent velocity, 
and its Rossby number can be as large as 10 given that its turbulent scale 
is about 0.3 kpc (see Poezd et al, 1993). Hence, the hot gas fills a quasi- 
spherical volume, where its pressure scale height of order 5 kpc is comparable 
to the disc radius. 

Ro = 1 at a scale 0.4 kpc in the warm gas, which is similar to the scale 
height of the gas layer. This implies that rotation significantly affects turbu- 
lent gas motions, making them helical on average. A convenient estimate of 
the associated a-effect can be obtained from Krause's formula, 

£ 2 fl 

a ~ — — 0.5 km s 1 , (3) 

where Q is the angular velocity, and the numerical estimate refers to the 
Solar neighbourhood of the Milky Way. Thus, ao ~ 0.05«o near the Sun and 
increases in the inner Galaxy together with Q. 

The spatial distribution of galactic rotation is known for thousands galax- 
ies (Sofue & Rubin, 2001) from systematic Doppler shifts of various spectral 
lines emitted by stars and gas. In this respect, galaxies are much better 
explored than any star or planet (including the Sun and the Earth) where 
reliable data on the angular velocity in the interior are much less detailed 
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Galactocentric distance, kpc Galactocentric distance, kpe 

Figure 3: (a): The rotation speed sQ(s) in the galactic midplane versus 
galactocentric radius s in the Milky Way (solid) (Clemens, 1985), and the 
generic Schmidt's rotation curve with Uq = 200 km s -1 , sq = 3 kpc and n = 1 
(dashed), (b): The corresponding rotation shear rates (taken with minus 
sign), — sdQ/ds. 



and reliable or even unavailable. The radial profile of the galactic rotational 
velocity is called the rotation curve. Rotation curves of most galaxies are flat 
beyond a certain distance from the axis, so Q oc s^ 1 is a good approximation 
for s >5kpc. The rotation curve of a generic galaxy, known as the Schmidt 
rotation curve and shown in Fig. El has the form 



sO(s) = U — 
■So 



1 + 2 [ _ 

3 T 3 1 



-i -3/2n 



where the parameters vary between various galaxies in the range so ~ 5- 
20 kpc, Uq ~ 200 km s _1 and n ~ 0.7-1. This rotation curve is not flat at large 
radii, but it provides an acceptable approximation at moderate distances 
from galactic centre where magnetic field generation is most intense. Some 
galaxies have more complicated rotation curves. Notably, the Milky Way and 
M31 are among them — see Figs 01 and |3J The complexity of the rotation 
curves is explained by a complicated distribution of the gravitating (stellar 
and dark) mass in those galaxies. It is evident from Fig. |3b that the rotation 
shear is strong at all radii even for the Schmidt rotation curve, and so the 
rotation in the inner part of a spiral galaxy cannot be approximated by the 
solid-body law, even if the shape of some rotating curves tempts to do so. 

The vertical variation of the rotation velocity is only poorly known. In 
a uniform gravitating disc of infinite radial extent the angular velocity of 
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rotation would be constant in z. Then it is natural to expect that Vt should 
decrease along z at a scale comparable to the radial scale length of the grav- 
itating mass in the disc, typically s* = 3-5 kpc. Recent observations of gas 
motions in galactic halos have confirmed such a decrease (Fraternali et ai, 
2003). In the absence of detailed models, an approximation Q oc exp (—z/s*) 
seems to be appropriate. 



3 Magnetic fields observed in galaxies 

Estimates of magnetic field strength in the diffuse interstellar medium of 
the Milky Way and other galaxies are most efficiently obtained from the 
intensity and Faraday rotation of synchrotron emission. Other methods are 
only sensitive to relatively strong magnetic fields that occur in dense clouds 
(Zeeman splitting) or are difficult to quantify (optical polarization of star 
light by dust grains). The total / and polarized P synchrotron intensities 
and the Faraday rotation measure RM are weighted integrals of magnetic 
field over the path length L from the source to the observer, so they provide 
a measure of the average magnetic field in the emitting or magneto-active 
volume: 

I = K n cr B 2 L ds , 

P = K J n^B^ds, (4) 
RM = Id j n c B\\ ds , 

where n cr and n e are the number densities of relativistic and thermal elec- 
trons, B is the total magnetic field comprising a regular B and random 
b parts, B = B + b with (B) = B, (b) = and (B 2 ) = B 2 + (b 2 ), 
angular brackets denote averaging, subscripts _L and || refer to magnetic 
field components perpendicular and parallel to the line of sight, and K and 
Ki = e 3 /(27rm^c 4 ) = 0.81 rad var 2 cm 3 /iG -1 pc _1 are certain dimensional 
constants (with e amd m e the electron charge and mass and c the speed of 
light). The degree of polarization p is related to the ratio (b 2 ) / B , 
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where the random field b has been assumed to be isotropic in the last equality, 
n CT is assumed to be a constant, and p ~ 0.75 weakly depends on the 
spectral index of the emission. This widely used relation is only approximate. 
In particular, it does not allow for any anisotropy of the random magnetic 
field, for the dependence of n CT on B, and for depolarization effects; some 
generalizations are discussed by Sokoloff et al. (1998). 

The orientation of the apparent large-scale magnetic field in the sky plane 
is given by the observed S-vector of the polarized synchrotron emission. Due 
to Faraday rotation, the true orientation can differ by an angle of RMA 2 , 
which amounts to 10°-20° at a wavelength A = 6 cm. The special importance 
of the Faraday rotation measure, RM, is that this observable is sensitive to 
the direction of B (the sign of Bn) and this allows one to determine not only 
the orientation of B but also its direction. Thus, analysis of Faraday rotation 
measures can reveal the three-dimensional structure of the magnetic vector 
field (Berkhuijsen et al, 1997; Beck et al., 1996). 

Since n CT is difficult to measure, it is often assumed that magnetic field and 
cosmic rays are in pressure equilibrium or energy equipartition; this allows 
to express n cr in terms of B. The physical basis of this assumption is the 
fact that cosmic rays (charged particles of relativistic energies) are confined 
by magnetic fields. An additional assumption involved is that the energy 
density of relativistic electrons responsible for synchrotron emission (energy 
of several GeV per particle) is one percent of the proton energy density in 
the same energy interval, as measured near the Earth. 

The cosmic ray number densityn cr in the Milky Way can be determined 
independently from 7-ray emission produced when cosmic ray particles inter- 
act with the interstellar gas. Then magnetic field strength can be obtained 
without assuming equipartition (Strong et al, 2000); the results are generally 
consistent with the equipartition values. However, Eq. (j3J) is not consistent 
with the equipartition or pressure balance between cosmic rays and magnetic 
fields as it assumes that n cr = const. Therefore, B obtained from Eq. (jSJ) can 
be inaccurate (Beck et al, 2003). 

The mean thermal electron density n e in the ISM can be obtained from 
the emission measure of the interstellar g observable defined as EM oc 
J L nl ds, but this involves the poorly known filling factor of interstellar clouds. 
In the Milky Way, the dispersion measures of pulsars, DM = J L n e ds provide 
information about the mean thermal electron density, but the accuracy is 
limited by our uncertain knowledge of distances to pulsars. Estimates of the 
strength of the regular magnetic field in the Milky Way are often obtained 
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Figure 4: The distributions of the strength of the total magnetic field in a 
sample of spiral galaxies obtained from the observed synchrotron intensity / 
using energy equipartition between magnetic fields and cosmic rays (p. 109 in 
Niklas, 1995) under slightly different assumptions. The estimates of the left- 
hand panel were derived from integrating the observed synchrotron intensity 
in the range corresponding to the relativistic electron energies from 300 MeV 
to infinity, and in the right-hand panel the integration was over a frequency 
range 10 MHz-10 GHz. Results presented in the left-hand panel are better 
justified physically (§2.1 in Beck et al, 1996; §III.A.l in Widrow, 2002). 



from the Faraday rotation measures of pulsars simply as 

This estimate is meaningful if magnetic field and thermal electron density are 
statistically uncorrelated. If the fluctuations in magnetic field and thermal 
electron density are correlated with each other, they will contribute positively 
to RM and Eq. © will yield overestimated B». In the case of anticorrelated 
fluctuations, their contribution is negative and Eq. © is an underestimate. 
As shown by Beck et al. (2003), physically reasonable assumptions about the 
statistical relation between magnetic field strength and electron density can 
lead to Eq. (jUJ) being in error by a factor of 2-3. 

The observable quantities (jlj) have provided extensive data on magnetic 
field strengths in both the Milky Way and external galaxies (Ruzmaikin et 
al.., 1988; Beck et al, 1996; Beck, 2000, 2001). The average total field 
strengths in nearby spiral galaxies obtained from total synchrotron inten- 
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sity I range from B « 4/iG in the galaxy M31 to about 15 /iG in M51, 
with the mean for the sample of 74 galaxies of B = 9 /iG (Beck, 2000). 
Figure 0] shows the distribution of magnetic field strength in a sample of 
spiral galaxies. The typical degree of polarization of synchrotron emission 
from galaxies at short radio wavelengths is p = 10-20%, so Eq. © gives 
B/B = 0.4-0.5; these are always lower limits due to the limited resolu- 
tion of the observations, and B/B = 0.6-0.7 is a more plausible estimate. 
Most existing polarization surveys of synchrotron emission from the Milky 
Way, having much better spatial resolution, suffer from Faraday depolar- 
ization effects and missing large-scale emission and cannot provide reliable 
values for p. The total equipartition magnetic field in the Solar neighbour- 
hood is estimated as B = 6 ± 2 p,G from the synchrotron intensity of the 
diffuse Galactic radio background (E. M. Berkhuijsen, in Beck, 2001). Com- 
bined with B/B = 0.65, this yields a strength of the local regular field of 
B — 4 ± 1 /iG. Hence, the typical strength of the local Galactic random 
magnetic fields, b = (B 2 — B 2 ) 1 / 2 = 5 ± 2 /iG, exceeds that of the regular 
field by a factor b/B = 1.3 ± 0.6. RM data yield similar values for this ratio 
(§IV.4 in Ruzmaikin et al, 1988). 

Meanwhile, the values of B in the Milky Way obtained from Faraday 
rotation measures seem to be systematically lower than the above values (see 
Beck et al, 2003, and references therein). RM of pulsars and extragalactic 
radio sources yield B = 1-2 /iG in the Solar vicinity, a value about twice 
smaller than that inferred from the synchrotron intensity and polarization. 
There can be several reasons for the discrepancy between the estimates of 
the regular magnetic field strength from Faraday rotation and synchrotron 
intensity. Both methods suffer from systematic errors due to our uncertain 
knowledge of thermal and relativistic electron densities, so one cannot be 
sure if the difference is significant. Nevertheless, the discrepancy seems to be 
worrying enough to consider carefully its possible reasons. 

The discrepancy can be explained, at least in part, if the methods de- 
scribed above sample different volumes. The observation depth of total syn- 
chrotron emission, starlight polarization and of Faraday rotation measures 
are all of the order of a few kpc. Polarized emission, however, may emerge 
from more nearby regions. However, a more fundamental reason for the dis- 
crepancy can be partial correlation between fluctuations in magnetic field 
and electron density. Such a correlation can arise from statistical pressure 
balance where regions with larger gas density have weaker magnetic field, and 



13 



vice versa. As discussed by Beck et al. (2003), the term (6||n e ) then differs 
from zero and contributes to the observed RM leading to underestimated 
B. In a similar manner, correlation between B and the cosmic ray num- 
ber density biases the estimates of magnetic field from synchrotron intensity 
and polarization (see also Sokoloff et al, 1998). Altogether, B = 4/zG and 
6 = 5 fiG seem to be acceptable estimates of magnetic field strengths near the 
Sun. The geometry and three-dimensional structure of the magnetic fields 
observed in spiral galaxies are further discussed in Sect. El 

4 The origin of galactic magnetic fields 

There are two basic approaches to the origin of global magnetic structures in 
spiral galaxies — one of them asserts that the observed structures represent a 
primordial magnetic field twisted by differential rotation, and the other that 
they are due to ongoing dynamo action within the galaxy. The simplicity of 
the former theory is appealing, but it fails to explain the strength, geome- 
try and apparent lifetime of galactic magnetic fields (Ruzmaikin et al., 1988; 
Beck et al, 1996; Kulsrud, 1999; Widrow, 2002; see Sect. El below). Fur- 
thermore, there are no mechanisms known to produce cosmological magnetic 
fields of required strength and scale (Beck et al, 1996), although Kulsrud et 
al. (1997) argue that suitable magnetic field can be produced in protogalaxies. 
Dynamo models appear to be much better consistent with the observational 
and theoretical knowledge of interstellar gas, and all models of magnetic fields 
in specific galaxies, known to the author, have been formulated in terms of 
dynamo theory. It seems to be very plausible that galactic magnetic fields 
are generated by some kind of dynamo action, i.e., that they are produced 
in situ. The most promising is the mean-field turbulent dynamo. 

4.1 Mean- field models of the galactic dynamo 

As discussed in Sect. 12.21 the discs of spiral galaxies are thin. This provides 
a natural small parameter, the disc aspect ratio, Eq. (J2J). This greatly facili- 
tates modelling of many global phenomena in galaxies, including large-scale 
magnetic fields. Parker (1971) and Vainshtein & Ruzmaikin (1971, 1972) 
were the first to suggest mean-field dynamo models for spiral galaxies. These 
were local models discussed in Sect. 14.1.21 where only derivatives across the 
disc (in z) are retained. The theory has been extended to two and more 
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dimensions and applied to specific galaxies (see Ruzmaikin et ai, 1988, Beck 
et ai, 1996, and Widrow, 2002, and references therein). Rigorous asymptotic 
solutions for the au;-dynamo in a thin disc were developed by Soward (1978, 
1992a,b) and further discussed by Priklonsky et al. (2000) and Willis et al. 
(2003). Reviews of these results can be found in Ruzmaikin et al. (1988), 
Beck et al. (1996), Kulsrud (1999) and Soward (2003). 

In this section we present asymptotic solutions of the mean-field dynamo 
equations in a thin disc surrounded by vacuum. We first consider axially sym- 
metric solutions of the kinematic problem, and then discuss generalizations 
to non-axisymmetric modes and to nonlinear regimes. Cylindrical coordi- 
nates (s, 0, z) with the origin the galactic centre and the z-axis parallel to 
the galactic angular velocity are used throughout this chapter. In this section 
we use dimensionless variables, with s and z measured in the units of the 
characteristic disc radius and disc half-thickness (e.g., sq = s & ~ 8.5 kpc and 
ho = 0.5 kpc), respectively. Then the dimensionless radial and axial distances 
are both of order unity within the disc as they are measured in different units 
in order to make the disc thinness explicit. The corresponding time unit is 
the turbulent magnetic diffusion time across the disc, h\jry\: ~ 7.5 x 10 s yr. 

It is convenient to introduce a unit rotational shear rate Go, 



with g(s,z) its dimensionless value and Go — — ^© for a flat rotation curve, 
Q oc s _1 , and adopt the characteristic magnitude of the a-coefficient near 
the Sun as given by Eq. (jSJ). 

4.1.1 Kinematic, axially symmetric solutions 

The three components of an axially symmetric magnetic field can be ex- 
pressed in terms of the azimuthal components of the large-scale magnetic 
field B,!, and vector potential A^, 



The dimensionless governing equations, resulting from standard mean-field 
dynamo equations have the form 



G = s— = G g(s,z) , 
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Vt Vt 
are the turbulent magnetic Reynolds numbers that characterize the inten- 
sity of induction effects due to differential rotation and the mean helicity of 
turbulence, respectively We have neglected the vertical shear dVt/dz which 
can easily be restored, and assumed for simplicity that Vt = const. A term 
containing a has been neglected in Eq. (J7J) for the sake of simplicity (but can 
easily be restored), so the equations are written in the acu-approximation. 

The kinematic, axially symmetric asymptotic solution in a thin disc has 
the form 




Q(e 



-1/3 




where F is the growth rate, (£>, ^4) represent the suitably normalized local 
solution (obtained for fixed s), and Q is the amplitude of the solution which 
can be identified with the field strength at a given radius. 



4.1.2 The local solution 

The local solution (with s fixed) arises in the lowest order in e. Its governing 
equations, obtained from Eqs ((Zj) and (jHJ) by putting e = 0, contain only 
derivatives with respect to z, with coefficients depend on s as a parameter 
(hence, the notation of the arguments of b and a with semicolon separating 
z and s): 

8 A d 2 B 

7 {S)B = -R u g( 8 )-£ + —. (10) 

1 (s)A = R a a{ S) z)B + ^, (11) 

were 7(5) is the local growth rate. The boundary conditions often applied 
at the disc surface z = ±h(s) correspond to vacuum outside the disc. For 
axisymmetric fields and to the lowest order in e they are (see below) 

dA 

B = and — = at z = ±h(s) . (12) 
oz 
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Since a is an odd function of z, kinematic modes have either even (quadru- 
pole) or odd (dipole) parity, with the following symmetry conditions at the 
disc midplane (see, e.g., Ruzmaikin et ai, 1988): 

OB 

= and ^4 = at z = (quadrupole) , (13) 

or 

dA 

B = and — = at z = (dipole) . (14) 

oz 

In order to clarify the nature of the dynamo modes in a thin disc, here 
we consider an approximate solution of Eqs (fTU|) and (fTTj) in the form of 
expansion in free-decay modes B n {z) and A n [z) obtained for R a = R w = 0: 

n _cPBn d 2 A n 



OZ 1 UZ A 

where 7„ (< 0) is the decay rate of the nth mode. For the boundary condi- 
tions (fT2|) and (JT3J) that select quadrupolar modes, the resulting orthonormal 
set of basis functions is given by 




V2 



cos 




^ v^sin 

72n = 72n+l = -VT 2 (n + |) 

The free-decay eigenvalues are all doubly degenerate, and two vector eigen- 
functions, one with odd index and the other with even one, correspond to 
each eigenvalue, one with &2n+i = 0, and the other with am = 0. The 
eigenfunctions are normalized to have Jq(B 2 + A\) dz = 1. 
The solution of Eqs (fTU|) and (fTTj) is represented as 




where c n are constants. We substitute this series into Eqs (fTUj) and (fTTj) . 
multiply by (bk, clu) and integrate over z from to h to obtain an algebraic 
system of homogeneous equations for whose solvability condition yields 
an algebraic equation for 7. For our current purposes, it is sufficient to 
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retain the smallest possible number of modes, which results in a system of 
two equations for cq and c\ and a quadratic equation for 7 whose positive 
solution is given by 



To assess the accuracy of Eq. (|15|). we note that it yields 7 = for D = 
D cr = — 7r 5 /16 ~ —19, as compared with the accurate value of D cr = —8 
(Ruzmaikin et al, 1988). This solution indicates that the dominant mode is 
non-oscillatory (Im7 = 0); this is confirmed by other analytical and numer- 
ical solutions of the dynamo equations in thin discs. 

A similar solution can be obtained for dipolar modes. The free decay 
modes of dipolar symmetry have 7 n = —n 2/ K 2 , n = 1,2,..., so that the 
lowest dipolar mode decays four times faster than the lowest quadrupolar 
mode. The reason for that is that the azimuthal field of dipolar parity has 
zero not only at \z\ = h but also at z = and so a smaller scale than the 
quadrupolar solution. This immediately implies that quadrupolar modes, 



with B^z) = B^-z), B s (z) = -B 8 (-z), B z {z) = B z (-z), should be 



dominant in galactic discs. The dominant symmetry of galactic magnetic 
fields is thus expected to be different from that in stars and planets, where 
dipolar fields are preferred. This prediction is confirmed by observations (see 
Sect. El). 

4.1.3 The global solution 

The vacuum boundary conditions are often used in analytical and semi- 
analytical studies of disc dynamos because of their (relative) simplicity. Most 
importantly, they have a local form in the lowest order in e — see Eq. (|12|1. 
However, this advantage is lost as soon as the next order in e is considered, 
which is needed in order to obtain a governing equation for the field distribu- 
tion along radius, Q. To this order, non-local magnetic connection between 
different radii has to be included, i.e., the fact that magnetic lines leave the 
disc at some radius, pass through the surrounding vacuum and return to the 
disc at another radius. In this section we discuss the radial dynamo equation, 




(15) 



where 




for a = simrz/h , 
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and for this purpose we have to consider vacuum boundary conditions to the 
first order in e. 

If the disc is surrounded by vacuum, there are no electric currents outside 
the disc, i.e., V x B = 0, so that the outer magnetic field is potential, B = 
— V$. Then axial symmetry implies that the azimuthal field vanishes outside 
the disc. Since magnetic field must be continuous on the disc boundary, this 
yields the following boundary condition at the disc surface z = ±h(s): 

B^\ z=±h = . (16) 

The vacuum boundary condition for the poloidal field (determined by A^) 
was derived in local Cartesian coordinates by Soward (1978). Priklonsky et 
al. (2000) rederived it in cylindrical geometry in the form 

^±- £ -C(A^)=0 atz = ±h(s), (17) 
where the integral operator C{A^) is defined as 

£ W=f (7^) ds ' 

with the kernel _ 

poo 

W(s, s') = ss' / J 1 (k8)J 1 (ks') dk , 
Jo 

where J\(x) is the Bessel function. Willis et al. (2003) obtained another, 
equivalent form of the integral operator involving Green's function of the 
Neumann problem for the Laplace equation. 

The integral part of the boundary condition (|1 7(1 can be transferred into 
a non-local term in the equation for Q which then becomes an integro- 
differential equation of the form (Priklonsky et al., 2000) 

[T- 1 (s)]q(s)=ep(s)C{q(s)} , (18) 

where 

q (s)=Q(s)A(h;s), K^ # ff S) " X = 

Here X is the eigenvector of the lowest-order boundary value problem dis- 
cussed in Sect. 14.1.21 the asterisk denotes the eigenvector of its adjoint prob- 
lem, and 

(X,X*> = / X-X*dz. 
Jo 
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The solution of Eq. (JTHJ) subject to the boundary conditions 

q(0) = and q — > as s — > oo 

provides yet another eigenvalue problem, for which the eigenvalue is the 
global growth rate T and the eigenfunction is q(s) which determines the 
radial profile of the global eigenfunction Q. As shown by Willis et al. (2003), 
the effect of the integral term in Eq. (fTHj) can be described as enhanced radial 
diffusion. 

Equation (|18|) is complicated enough as to provoke an irresistible desire 
to simplify it. Such a simplification, employed by Baryshnikova et al. (1987) 
(see also Ruzmaikin et al., 1988) consists of neglecting the term containing 
e in the boundary condition ([17)1 . This makes the boundary condition local 
and leads to the following equation for Q(s): 

' r = 4 (19) 

similar to Eq. 1)18)1. but with the integral term replaced by the diffusion 
operator. Formally, Eq. (fT9"|) can be obtained from Eq. (|T%j) by replacing the 
integral kernel by the delta-function, W(s,s') — > 5(s — s'). In other words, 
this simplification neglects any nonlocal coupling between different parts of 
the disc via the halo, but includes the local diffusive coupling within the 
disc. We note in this connection that the kernel W(s, s') is indeed singular, 
although the singularity is only logarithmic in reality, W(s, s') ~ In \s — s'\. 

The above simplification greatly facilitates the analysis of the global dy- 
namo solutions and all applications of the thin-disc asymptotics to galaxies 
and accretion discs neglect the nonlocal effects. Equation (fHJj) can be readily 
solved using a variety of analytical and numerical techniques (Ruzmaikin et 
al., 1988), but some features of the solution are lost together with nonlocal 
effects. The most important failure is that the asymptotic scaling of the 
solution with e is affected, with the radial scale becoming e -1//2 /io instead of 
the correct value £ _1 / 3 /io- However, the difference is hardly significant nu- 
merically for the realistic values e ~ 10 -1 -10~ 2 . We note that the thin-disc 
asymptotics are reasonably accurate for e < 10 _1 (Baryshnikova et al, 1987; 
Willis et al, 2003). 

Another effect of the nonlocal effects is that solutions of Eq. (JTHj) possess 
algebraic tails far away from the dynamo active region, q ~ s -4 , whereas 
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solutions of Eq. (J 19)) have exponential tails typical of the diffusion equa- 
tion. This affects the speed of propagation of magnetic fronts during the 
kinematic growth of the magnetic field: with the nonlocal effects, the fronts 
propagate exponentially, whereas the local radial diffusion alone results in a 
linear propagation. 

These topics are discussed in detail by Willis et al. (2003) who compare 
numerical solutions of Eqs (|18|) and (|19j) . Whether or not the nonlocal ef- 
fects can be neglected depends on the goals of the analysis. There are several 
reasons why this simplification appears to be justified. The neglect of non- 
local effects does not seem to affect significantly any observable quantities, 
whereas the parameters of spiral galaxies and of their magnetic fields are 
known with a rather limited accuracy anyway. Moreover, the halos of spiral 
galaxies can be described as vacuum only in a very approximate sense, and 
the finite conductivity of the halo will weaken the nonlocal effects. 

4.1.4 Non-axisymmetric, nonlinear and numerical solutions 

The above asymptotic theory can readily be extended to non-axisymmetric 
solutions. This generalization is discussed by Krasheninnikova et al. (1989) 
and Ruzmaikin et al. (1988). Starchenko & Shukurov (1989) developed WKB 
asymptotic solutions of the mean-field galactic dynamo equations valid for 
\D\ ^> 1. A similar asymptotic regime for one-dimensional dynamo equations 
Ifflty and (HD is discussed in §9.IV of Zeldovich et al. (1983). 

Another useful approximate approach, known as the 'no-z' approxima- 
tion, was suggested by Subramanian & Mestel (1993). In this approxima- 
tion, derivatives across the disc in Eqs. (J7J) and (jHI) or their three-dimensional 
analogues are replaced by division by the disc semi-thickness, d/dz —>■ 1/h, 
and the resulting equations in s and <fi are solved, e.g., by the WKB method 
or numerically. This approach appears to be rather crude at first sight, but 
it is quite efficient because the structure of the magnetic field across a thin 
disc is quite simple, at least for the lowest mode. A refinement of the ap- 
proximation to improve its accuracy is discussed by Phillips (2001). Mestel 
and Subramanian (1991) and Subramanian & Mestel (1993) apply these so- 
lutions to study the effects of spiral arms on galactic magnetic fields. This 
approximation was also extensively used in numerical simulations of galactic 
dynamos (Moss 1995; see Moss et al, 2001 for an example). 

Nonlinear asymptotics of Eqs (fTU|) and (|TT|) for \D\ 3> 1 are discussed 
by Kvasz et al. (1992), where it is supposed that the nonlinearity affects 
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significantly magnetic field distribution across the disc, and to the lowest 
approximation the steady state of the dynamo is established locally. This, 
however, may not be the case. The radial coupling is significant already at the 
kinematic stage where it results in the establishment of a global eigenfunction 
as described by Eq. (j!8|) or Eq. (fT9"|) . Nonlinear effects are more likely to affect 
the global eigenfunction, and so have to affect the radial equation. Poezd et 
al. (1993) have derived a nonlinear version of Eq. (J19|) assuming the standard 
form of a-quenching with the a-coefficient modified by magnetic field as 

a = — 2 , 20 

1 + B 2 /Bl 

where Bq is a suitably chosen saturation level most often identified with a 
state where magnetic and turbulent kinetic energy densities are of the same 
order of magnitude. As a result, the magnetic field can grow when B <C B , 
but then the growth slows down as the quenched dynamo number obtained 
with a approaches its critical value D cr , and the field growth saturates at 
B pa Bq. In terms of the thin-disc asymptotic model, this implies that 7(5) 
in Eqs ffTS|) and (fTUjl ought to be replaced by 7(5) (1 — Q 2 /Bq), so that the 
nonlinear version of Eq. (J 19)) with the nonlinearity (|2U|) has been derived in 
the form 

d 



d ® ( W 1 




S OS 



(21) 



provided the local solution has been normalized in such a way that Q is a 
field strength averaged across the disc at a given radius. The derivation of 
this equation by averaging the governing equations across the disc can be 
found in Poezd et al. (1993). This equation and its nonaxisymmetric version 
have been extensively applied to galactic dynamos (see Beck et al, 1996, and 
references therein). 

The detailed physical mechanism of the saturation of the dynamo ac- 
tion is still unclear. Cattaneo et al. (1996) suggest that the saturation is 
associated with the suppression of the Lagrangian chaos of the gas flow by 
the magnetic field. This mechanism, attractive in the context of convective 
systems (where the flow becomes random due to intrinsic reasons, e.g., insta- 
bilities), can hardly be effective in galaxies where the flow is random because 
of the randomness of its driving force (the supernova explosions). 

Most numerical solutions of galactic dynamo equations that extend be- 
yond the thin-disc approximation rely on the 'embedded disc' approach (Ste- 
pinski & Levy, 1988; Elstner et al, 1990). Instead of using complicated 
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boundary conditions at the disc surface, this approach considers a disc em- 
bedded into a halo whose size is large enough as to make unimportant bound- 
ary conditions posed at the remote halo boundary Since turbulent magnetic 
diffusivity in galactic halos is larger than in the disc (Sokoloff & Shukurov, 
1990; Poezd et al, 1993), meaningful embedded disc models are compatible 
with thin-disc asymptotic solutions obtained with vacuum boundary condi- 
tions and confirm the asymptotic results. The embedded disc approach was 
also used to study dynamo-active galactic halos (Brandenburg et al, 1992, 
1993, 1995; Elstner et al, 1995). Further extensions of disc dynamo models 
include the effects of magnetic buoyancy (Moss et al, 1999), accretion flows 
(Moss et al, 2000) and external magnetic fields (Moss & Shukurov, 2001, 
2004). 

An implication of the nonlinear model for the thin-disc dynamo is that 
the local solution is unaffected by nonlinear effects whose main role is to 
modify the radial field structure. An important consequence of this is that 
it can be reasonably expected that the pitch angle of magnetic lines, p = 
aictanBg/B^, is weakly affected by nonlinear effects, and so represents an 
important feature of the solution that can be directly compared with obser- 
vations (Baryshnikova et al, 1987). This expectation seems to be confirmed 
by observations fSect. IB~T|) . Nevertheless, the modification of the magnetic 
pitch angle by nonlinear effects has never been studied in detail, which seems 
to be a regrettable omission. 

4.1.5 Dynamo control parameters in spiral galaxies 

A remarkable feature of spiral galaxies is that they are (almost) transparent 
to electromagnetic waves over a broad range of frequencies, so the kinematics 
of the ISM is rather well understood, and therefore all parameters essential 
for dynamo action are well restricted by observations. This leaves less room 
for doubt and less freedom for speculation than in the case of other natural 
dynamos. Another advantage is that observations of polarized radio emission 
at a linear resolution of 1-3 kpc (typical of the modern observation of nearby 
galaxies) reveal exactly that field which is modelled by the mean-field dynamo 
theory (given volume and ensemble averages are identical). 

The mean-field dynamo is controlled by two dimensionless parameters 
quantifying the differential rotation and the so-called a-effect, as defined in 
Eq. @. Using Eqs. and Q and assuming a flat rotation curve, Q — Uq/s, 



23 



we obtain the following estimates for the Solar vicinity of the Milky Way: 

Uq IS Uq S 

where Uq = sqQq is the typical rotational velocity Since \Ru\ 3> R a , differen- 
tial rotation dominates in the production of the azimuthal magnetic field (i.e., 
the au;-dynamo approximation is well applicable), and the dynamo action is 
essentially controlled by a single parameter, the dynamo number 

D = R a R„ ~ 1(4 sQ ^ « -10 ( 2 « -20 , (22) 

Uq OS \ UqS J 

where the numerical estimate refers to the Solar vicinity. Thus, \D\ does 
exceed the critical value for the lowest, non-oscillatory quadrupole dynamo 
mode, which then can be expected to dominate in the main parts of spiral 
galaxies. It is often useful to consider the local dynamo number D(s), a func- 
tion of galactocentric radius s, obtained when the s-dependent, local values 
of the relevant parameters are used in Eq.|J3 instead of the characteristic ones. 

The local regeneration (e-folding) rate of the regular magnetic field 7 is 
related to the magnetic diffusion time along the smallest dimension of the 
gas layer and to the dynamo number (if \R^\ 3> R a ). Using the perturbation 
solution of Sect. 14.1^ the following expression (written in dimensional form) 
can be used as a rough estimate: 

7~^(y/\D\-y/\D^, tov\D\>D CT , (23) 

where D CT « —6 can be adopted from the more accurate numerical solution 
and numerical factor of order unity has been omitted. This yields the local 
e-folding time 7 _1 ~ 4 x 10 8 yr for the Solar neighbourhood. When the 
radial diffusion is included, i.e., Eq. (|TH|) is solved, the growth rate decreases, 
yielding a global e-folding time of ~ 10 9 yr near the Sun. Thus, the 
large-scale magnetic field near the Sun can be amplified by a factor of about 
10 4 during the galactic lifetime, 10 10 yr, and the Galactic seed field had to 
be rather strong, about 10~ 10 G. The fluctuation dynamo can produce such 
a statistical residual magnetic field at the scale of the leading eigenfunction 
either in the young galaxy (§VTI.13 in Ruzmaikin et ai, 1988; Widrow, 2002) 
or in the protogalaxy (Kulsrud et ai, 1997). 
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The above growth rate, estimated for the Solar neighbourhood of the 
Milky Way, is often erroneously adopted as a value typical of spiral galaxies in 
general. It is then important to note that the regeneration rate is significantly 
larger in the inner Galaxy (the local dynamo number rapidly grows as s 
becomes smaller, D(s) oc GQ - - see Fig. and in other galaxies. For 
example, Baryshnikova et al. (1987) estimate the global growth time of the 
leading axisymmetric mode in the galaxy M51 as 5 x 10 7 yr. 

Gaseous discs of spiral galaxies are flared, ie., h oc s + const at s > lOkpc, 
whereas uq only slightly varies with s. For a flat rotation curve, Q oc 
Eq. (|22JI then shows that the local dynamo number does not vary much with 
galactocentric radius s and remains supercritical, |-D(s)| > |-D cr | out to a 
large radius. It is therefore not surprising that regular magnetic fields have 
been detected in all galaxies where observations have sufficient sensitivity 
and resolution (Wielebinski & Krause, 1993; Beck et al, 1996; Beck, 2000, 
2001). 

A standard estimate of the steady-state strength of magnetic field pro- 
duced by the mean-field dynamo follows from the balance of the Lorentz 
force due to the large-scale magnetic field and the Coriolis force that causes 
deviations from mirror symmetry (Ruzmaikin et al, 1988; Shukurov, 1998): 

D \1 1/2 

1) (24) 
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where p ~ 1.7 x 10~ 24 gcm -3 is the density of interstellar gas and n its 
number density, n = p/mu with mn the proton mass. This estimate yields 
values that are in good agreement with observations, but its applicability 
has to be reconsidered in view of the current controversy about the nonlinear 
behaviour of mean-field dynamos (see Sect. 14.3)) . 

It is now clear what information is needed to construct a useful dynamo 
model for a specific galaxy: its rotation curve, the scale height of the gas layer, 
the turbulent scale and speed, and the gas density. All these parameters are 
observable, even though their observational estimates may be incomplete 
or controversial. One of successes of the mean-field dynamo theory is its 
application to spiral galaxies, where even simplest, quasi-kinematic models 
presented above are able to reproduce all salient features of the observed 
fields, both in terms of generic properties and for specific galaxies (Ruzmaikin 
et al, 1988). We discuss this in Sect.EJ 
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Recent observational progress has allowed to explore the effects of galactic 
spiral patterns on magnetic fields (Beck, 2000). The corresponding dynamo 
models require the knowledge of the arm-interarm contrast in all the relevant 
variables (Shukurov & Sokoloff, 1998; Shukurov, 1998; Shukurov et al, 2004). 

4.2 The fluctuation dynamo and small-scale magnetic 
fields 

Similarly to mean-field dynamos, the theory of the fluctuation dynamo is 
well understood in the kinematic regime, but nonlinear effects remain con- 
troversial. In this section we present results obtained with kinematic models 
of the fluctuation dynamo and those derived with simplified nonlinearity. 
The pioneering kinematic model of the fluctuation dynamo was developed 
by Kazantsev (1967), and many more recent developments are based on it. 
Detailed reviews of the theory and references can be found in §8. IV of Zel- 
dovich et al. (1983), Ch. 9 of Zeldovich et al. (1990) and in Brandenburg & 
Subramanian (2004). 

The growth time of the random magnetic field in a random velocity field 
of a scale £ is as short as the eddy turnover time, £/uq ~ 10 7 yr in the warm 
phase for £ = 0.1 kpc. The magnetic field produced by the dynamo action 
is a statistical ensemble of magnetic flux ropes whose length is of the order 
of the flow correlation length, £ ~ 0.05-0.1 kpc. The rope thickness is of the 
order of the resistive scale, £Rm _1//2 , in a single-scale velocity field, where Rm 
is the magnetic Reynolds number. A phenomenological model of dynamo in 
Kolmogorov turbulence yields the rope thickness of £Rm~ 3//4 (Subramanian, 
1998). The dynamo action can occur provided Rm > Rm cr , where the critical 
magnetic Reynolds number is estimated as Rm cr = 30-100 in simplified mod- 
els of homogeneous, incompressible turbulence. Recent studies have revealed 
the possibility that small-scale magnetic fields can have peculiar fine struc- 
ture because the magnetic dissipation scale in the interstellar gas is much 
smaller than that of turbulent motions, i.e., because the magnetic Prandtl 
number is much larger than unity (Schekochikhin et al., 2002). 

Subramanian (1999) suggested that a steady state, reached via the back- 
action of the magnetic field on the flow, can be established by the reduction 
of the effective magnetic Reynolds number down to the value critical for the 
dynamo action, an idea similar to the concept of a-quenching in the mean- 
field theory. Then the thickness of the ropes in the steady state can be 
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estimated as ffim" ' or £Rm~ ' . Using a model nonlinearity in the induc- 
tion equation with incompressible velocity field, Subramanian (1999) showed 
that the magnetic field strength within the ropes bo saturates at the equipar- 
tition level with kinetic energy density, &o/87r ~ \pui- The average magnetic 
energy density is estimated as b 2 /8n ~ ^Rm" 1 pu^, implying the volume fill- 
ing factor of the ropes of order f v ~ Rm^ 1 ~ 0.01. Correspondingly, the 
mean magnetic energy generated by the small-scale dynamo in the steady 
state is about 1% of the turbulent kinetic energy density, in agreement with 
numerical simulations. 

Shukurov & Berkhuijsen (2003) interpret thin, random filaments of zero 
polarized intensity observed in polarization maps of the Milky Way (known 
as canals) as a result of Faraday depolarization in the turbulent interstellar 
gas. This interpretation has resulted in a tentative estimate of the Taylor 
microscale of the interstellar turbulence 

-1/2 

£ T = £Rm ~0.6pc, 

where Rm is the effective magnetic Reynolds number in the ISM. This yields 
the following estimate: 

Rm =~ 10 4 . 

Of course, this is a very tentative estimate, and further analyses of observa- 
tions and theoretical developments will be needed to refine it. The value of 
Rm obtained is significantly larger than Rm cr obtained in idealized models. 
This might be due to the transonic nature of interstellar turbulence as the 
gas compressibility appears to hinder dynamo action. Kazantsev et al. (1985) 
have shown that the e-folding time of magnetic field in the acoustic-wave tur- 
bulence (i.e., a compressible flow) is as long as M^l/uo, where M. ( > 1) is 
the Mach number. 

Using parameters typical of the warm phase of the ISM, this theory pre- 
dicts that the small-scale dynamo would produce magnetic flux ropes of the 
length^ or the curvature radius) of about i = 50-100 pc and thickness 5-10 pc 
for Rm = 10 2 and 0.5-10 pc for Rm = 10 4 . The field strength within the 
ropes, if at equipartition with the turbulent energy, has to be of order 1.5 /zG 
in the warm phase (n = 0.1 cm -3 , uq = lOkrns -1 ) and 0.5 fiG in the hot gas 
(n = 10 _3 cm~ 3 , Uq = 40kms _1 ). Note that some heuristic models of the 
small-scale dynamo admit solutions with magnetic field strength within the 
ropes being significantly above the equipartition level, e.g., because the field 
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configuration locally approaches a force-free one, |(V x B) x B| <C B 2 /£, 
where i is the field scale (Belyanin et ai, 1993). 

The small-scale dynamo is not the only mechanism producing random 
magnetic fields (e.g., §4.1 in Beck et ai, 1996, and references therein). Any 
mean-field dynamo action producing magnetic fields at scales exceeding the 
turbulent scale also generates small-scale magnetic fields. Similarly to the 
mean magnetic field, this component of the turbulent field presumably has 
a filling factor close to unity in the warm gas and its strength is expected 
to be close to equipartition with the turbulent energy at all scales. This 
component of the turbulent magnetic field may be confined to the warm gas, 
the site of the mean-field dynamo action, so magnetic field in the hot phase 
may have a better pronounced ropy structure. 

The overall structure of the interstellar turbulent magnetic field in the 
warm gas can be envisaged as a quasi-uniform fluctuating background with 
one percent of the volume occupied by flux ropes (filaments) of a length 50- 
100 pc containing a well-ordered magnetic field. This basic distribution would 
be further complicated by compressibility, shock waves, MHD instabilities 
(such as Parker instability), the fine structure at subviscous scales, etc. 

The site of the mean-field dynamo action is plausibly the warm phase 
rather than the other phases of the ISM. The warm gas has a large filling 
factor (so it can occupy a percolating global region), it is, on average, in a 
state of hydrostatic equilibrium, so it is an ideal site for both the small-scale 
and mean-field dynamo action. Molecular clouds and dense Hi clouds have 
too small a filling factor to be of global importance. Fletcher & Shukurov 
(2001) argue that, globally, molecular clouds can be only weakly coupled 
to the magnetic field in the diffuse gas, but Beck (1991) suggests that a 
significant part of the large-scale magnetic flux can be anchored in molecular 
clouds. The time scale of the small-scale dynamo in the hot phase is £/uq ~ 
10 6 yr for uq = 40kms _1 and I = 0.04 kpc (the width of the hot, 'chimneys' 
extended vertically in the disc). This can be shorter than the advection 
time due to the vertical streaming, h/U z ~ 10 7 yr with h = lkpc and U z = 
100 km s -1 . Therefore, the small-scale dynamo action should be possible in 
the hot gas. However, the growth time of the mean magnetic field must be 
significantly longer than £/uq, reaching a few hundred Myr. Thus, the hot 
gas can hardly contribute significantly to the mean-field dynamo action in 
the disc and can drive the dynamo only in the halo (Sokoloff & Shukurov, 
1990). The main role of the fountain flow in the disc dynamo is to enhance 
magnetic connection between the disc and the halo (see Sect. 14.3)) . 



28 



4.3 Magnetic helicity balance in the galactic disc 

Conservation of magnetic helicity \= (A ■ B) (where B = V x A) in a per- 
fectly conducting medium has been identified as an important constraint on 
mean-field dynamos that plausibly explains the catastrophic quenching of the 
a-effect discussed elsewhere in this volume (Blackman & Field, 2000; Kleeorin 
et al, 2000, 2003; Brandenburg & Subramanian, 2004). In a closed system, 
magnetic helicity can only evolve on the (very long) molecular diffusion time 
scale; in galaxies, this time scale by far exceeds the Hubble time. The large- 
scale galactic magnetic fields have significant magnetic helicity of the order 
of LBgB^ ~ — ^LB 2 , where L > lkpc is the field scale, Bg/B^ = ship with 
p pa —15° the magnetic pitch angle. Since the initial (seed) magnetic field 
was weak, and so had negligible magnetic helicity, the large-scale magnetic 
helicity in a closed system must be balanced by the small-scale helicity of the 
opposite sign, ~ £hb 2 , where t\ v is an appropriate dominant scale of magnetic 
helicity. This immediately results in an upper limit on the steady-state mean 
magnetic field (Brandenburg & Subramanian, 2004, and references therein) 

= <4^0.4, (25) 

where the numerical value is obtained for £^ = 0.1 kpc and L = 1 kpc. The 
result of Vainshtein & Cattaneo (1992), B /b 2 ~ Rm _1//2 is recovered for 
£h ~ LRm" 1 / 2 . The observed relative strength of the mean field in spiral 
galaxies is given by B jb 2 ~ 0.5. The upper limit on the strength of the mean 
magnetic field (J2*5j) appears to be much lower than the observed field only 
if £h <C 0.1 kpc. For £ h = 0.1 kpc, the observed field strength is compatible 
with magnetic helicity conservation. What is perhaps more worrying, is that 
the mean magnetic fields can only grow at the long molecular diffusion time 
scale to reach this strength. 

Blackman & Field (2000) and Kleeorin et al. (2000) suggested that the 
losses of the small-scale magnetic helicity through the boundaries of the dy- 
namo region play the key role in the mean-field dynamo action. This is an 
appealing idea, especially because the mean-field dynamos rely on magnetic 
flux loss through the boundaries (§9.11 in Zeldovich et al, 1983; §VII.5 in 
Ruzmaikin et al, 1988). A similar situation occurs with the magnetic mo- 
ment, which is a conserved quantity, and it only grows in a dynamo system of 
a finite size because the dynamo just redistributes it expelling magnetic mo- 
ment out from the dynamo active region (Moffatt, 1978). However, these are 
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the mean magnetic flux and moment that need to be transferred through the 
boundaries. Transport by turbulent magnetic diffusion is sufficient for these 
purposes. The new aspect of the magnetic helicity balance is that healthy 
mean-field dynamo action requires asymmetry between the transports of the 
magnetic helicities of the large- and small-scale magnetic fields. 

A useful framework to assess the effects of magnetic helicity flow through 
the boundaries of the dynamo region was proposed by Brandenburg et al. 
(2002) who have presented the balance equation of magnetic helicity in the 
form 

^f + ^f = ~ 2 VXJ ~ ~Qb ~Qb , (26) 

where xb = A-B and Xb = ( a 'b) are the magnetic helicities of the mean and 
random magnetic fields, respectively, 77 is the molecular magnetic diffusivity, 
Xj = J ■ B and Xj = (j ' b) are the current helicities (with J = V x B 
the current density) . The first two terms on the right-hand side of Eq. (J26|) 
are responsible for the Ohmic losses whereas the last two terms represent 
the boundary losses. For illustrative purposes and following Brandenburg et 
al. (2002), we adopt the following assumptions, (i) The magnetic fields are 
fully helical, so M B = \k B \\B\ and M b = ^kb\xb\, where M B and M b are 
the average energy densities of the mean and random magnetic fields and k B 
and kb are their wave numbers, respectively. Furthermore, xj = k B XB and 
Xj = k\xb- (h) The mean and random magnetic fields have widely separated 
scales, k B -C kb- (iii) Approximate equipartition is maintained between the 
mean and random magnetic fields, M B « Mb. Then 



Xb 



Xb 



k b M B = kj xj 
k B M b k% Xj 



and so \xb\ 3> \Xb\ and \xj\ <C \Xj\- Assuming for definiteness that xbi Xj > 
0, we have Xb, Xj < 0, and Eq. can be approximated by 



dM B 
dt 



^k h k B M h + \k B {Q B + Q b ) . (27) 



It is important to note that the effective advection velocities for the large- 
scale and small-scale magnetic fields are not equal to each other. Both small- 
scale and large-scale magnetic fields are advected from the disc by the galactic 
fountain flow. With a typical vertical velocity of order U z = 100-200 km s -1 
and the surface covering factor of the hot gas / = 0.2-0.3, the effective 
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vertical advection speed is fU z = 30-70 km s -1 . However, the large-scale 
magnetic field is subject to turbulent pumping (turbulent diamagnetism). 
Given that the turbulent magnetic diffusivity in the disc and the halo are 
i]^ = 10 26 cm 2 s _1 and 77^ = 2 x 10 27 cm 2 s _1 (Poezd et al, 1993), respec- 
tively, and that the transition layer between the disc and the halo has a thick- 
ness of Az = lkpc, the resulting advection speed is U d = — |V?7t ~ — (2- 
3) km s _1 . Thus, the vertical advection velocities of the large-scale and small- 
scale magnetic fields are fU z + Ud and fU z , respectively. 

Now we can estimate the magnetic helicity fluxes through the disc surface 

as 

Qb — — (u B + j—] M B , Q b = U b M b , 

where = l/(4r]Tk B ) is the time scale of the (turbulent) diffusive transport 

of the mean magnetic field through the boundary, and U b and U b are effective 

advection velocities for the large-scale and small-scale magnetic helicities, 

respectively. The latter can be estimated from the following arguments. 

Consider advection of magnetic field through the disc surface z = h by a flow 
2 2 

with a speed U, OB jdt = —UdB jdz. Assuming for simplicity that U is 

2 

independent of z, we obtain by integration over z: 2hM B = —2U B (h), where 

M B = (2h)~ l j\~B dz. With M B = -k B x B /2, this shows that advection of 

magnetic field at a speed U produces the large-scale helicity loss at a rate 

Xb = Q B = (2U/kBh)B (h). Here B(h) is the large-scale field strength at 

2 

the disc surface, which is given by B (h) = (,M B , where ^ < 1 because the 
large-scale magnetic field at the surface must be weaker than that deep in the 
disc. For example, £ 1 for vacuum boundary conditions where B^h) = 0. 
Thus, 

U B = ^r(fU z + U d ) . 
k B h 

Unlike the large-scale magnetic field, the small-scale magnetic fields are not 
necessarily weaker at the disc surface, so similar arguments yield 

U b = ^-fU z ± U B • 
k b h 

Thus, there are several reasons for the magnetic helicity fluxes through 
the disc surface to be different at small and large scales: most importantly, 
the large-scale magnetic field at the surface can be much smaller than that 
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deep in the disc (£ <C 1) and, in addition, turbulent diamagnetism introduces 
further difference (U& ^ 0). 

Equation (j2*Tj) has the following solution satisfying the initial condition 
M B (0) = 0: 



M B A V k b + U b 



Ar] T k B + U B 



exp 



1 



k B U B \ t 



(28) 



For t <C tb, this solution captures the exponential growth of the mean 
magnetic field at a time scale 2tb, Mb oc t/2r B - For U B = U b = 0, 
we obtain Ms/M b m r]k b /r)TkB ~ Rrn -1 — this corresponds to the catas- 
trophic quenching of the a-effect associated with approximate magnetic he- 
licity conservation in a medium with (weak) Ohmic losses alone. However, 
for U b S> 4:Tjk b (a condition safely satisfied for any realistically small if) and 
Ub 3> 477t&b ~ 8kms -1 , we obtain 



M f 



M h 



(29) 



t— »oo 



Ub_ ks_ ±_ 
U B ~ ^6 ~ 10^ ' 

where we recall that £ < 1 and neglect L^. Thus, states with M B w M 6 
cannot be excluded, and this equipartition state is reached at the time scale 
of order tb ~ 4 x 10 8 yr. 

These arguments suggest that the growth rate of the mean magnetic field 
is limited from above by the flux of the mean magnetic helicity through the 
boundary of the dynamo region, whereas the upper limit for its steady state 
strength is controlled by the rate at which the small-scale magnetic helicity 
is transferred through the boundaries, Eq. ([29)1 . Another limit on the mean 
field strength arises from the balance of the Lorentz and Coriolis forces in the 
disc, Eq. (I24J) . The steady-state strength of the mean magnetic field is the 
minimum of the two values. These arguments suggest that the restrictions 
on the mean-field dynamo action from magnetic helicity conservation can be 
removed as soon as one allows for the disc-halo connection and fountain flows 
in spiral galaxies. Of course, these heuristic arguments have to be confirmed 
by quantitative analysis. 
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Figure 5: Left panel: The pitch angle of magnetic field in the galaxy M31 
as obtained from radio polarization observations (squares with error bars) 
(Fletcher et al, 2004), from Eq. (|3*Tj) using the rotation curve of Deharveng & 
Pellet (1975) and Haud (1981) (dashed) and Braun (1991) (dotted), and from 
Eq. (|3*3^1 with D cr = 1 using the same rotation curves (open and filled circles, 
respectively); h(r) is twice the Hi scale height of Braun (1991). Results from 
a nonlinear dynamo model for M31 (Beck et al., 1998) are shown with solid 
line. Right panel: The rotation curve of M31 from Deharveng Sz Pellet (1975) 
and Haud (1981) (solid) and from Braun (1991) (dashed). 
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5 Observational evidence for the origin of 
galactic magnetic fields 

5.1 Magnetic pitch angle 

Regular magnetic fields observed in spiral galaxies have field lines in the form 
of a spiral with a pitch angle in the range p = — (10°-30°), with negative 
values indicating a trailing spiral (e.g., Beck et ai, 1996). As discussed in 
Sect. 14.1741 the value of the pitch angle is a useful diagnostic of the mechanism 
maintaining the magnetic field. 

Consider the simplest from of mean-field dynamo equations (jlUj) and (fTTj) 
appropriate for a thin galactic disc, but now written in terms of dimensional 
variables for B s and Bf. 

dB s d d 2 B s dBj, d 2 B^ 

— -^ (a ^ )+r > T ^' ~w = GBs+ ^^- (30) 

Any regular magnetic field maintained by the dynamo must have a non-zero 
pitch angle: for B s = (a purely azimuthal magnetic field), equation for 
in (jHUj) reduces to a diffusion equation dB^/dt = ri^d 2 B^jdz 2 which only 
has decaying solutions, oc exp(— r] T t/h 2 ). The same applies to a purely 
radial magnetic field. 

Consider exponentially growing solutions, B s ^ = £><, ^exp(7t), and replace 
d/dz by 1/h and d 2 /dz 2 by —1/h 2 (as in the 'no- z' approximation) to obtain 
from Eqs. (J3Uj) two algebraic equations, 

(7 + m/h 2 ) B s + aB^/h = , -GB S + (7 + m/h 2 ) ^ = , 

which have non-trivial solutions only if the determinant vanishes, which yields 
(7 + VT/h 2 ) 2 ~ -aG/h, and Eq. follows with D cr = 1. The resulting 
estimate of the magnetic pitch angle is given by 



B s 

t&np = — 





R a 


£ 

- rsj 

h 




l-Gh~ \j 




dins 



-1/2 

(31) 



For ijh ~ 1/4 and a flat rotation curve, <91nf2/<91ns = —1, we obtain 
p ~ —15°, and this is the middle of the range observed in spiral galaxies. 
More elaborate treatments discussed by Ruzmaikin et al. (1988b) confirm 
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this estimate of p and yield a more accurate value of D CT . For example, the 
perturbation solution of Sect. 14. l."2l yields 



If the steady state is established by reducing R a to its critical value as 
to obtain R a Ru> = D CT , then the pitch angle in the nonlinear steady state 
becomes 



The magnetic pitch angle in M31 determined from observations and dy- 
namo theory is shown in Fig. Although the model curves show noticeable 
differences from the observed pitch angles, the general agreement is encour- 
aging. The situation is typical: magnetic pitch angles of spiral galaxies are 
in a good agreement with predictions of dynamo theory (Beck et ai, 1996). 

This picture does not explain why the pitch angles of galactic magnetic 
fields are invariably close (though not equal) to those of the spiral pattern 
in the parent galaxy. A plausible explanation is that magnetic pitch angles 
are further affected by streaming motions associated with the spiral pattern 
to make the match almost perfect (Moss, 1998). We note, however, that the 
pitch angle of the large-scale magnetic field near the Sun differs significantly 
from that of the local (Orion) arm; it is not clear whether this misalignment 
is of a local or global nature. 

As shown by Moss et al. (2000), magnetic pitch angle can be affected by 
an axisymmetric radial inflow (as well as outflow): 



which is useful to compare with Eqs (pTTj) and ((H2J)- This effect is important 
if u s >2r] T /h ~ lkms" 1 (cf. Sect.Q. 

Twisting of a horizontal primordial magnetic field by galactic differential 
rotation leads to a tightly wound magnetic structure with magnetic field 
direction alternating with radius at a progressively smaller scale As ~ s /\G\t 
with tanp ~ — where So ~ lOkpc is the scale of variation in Q 

(see §3.3 in Moffatt, 1978; Kulsrud, 1999 for a detailed discussion). The 
winding-up proceeds until a time t ~ 5 x 10 9 yr such that \G\t ~ IC^I 1 / 2 , 




(32) 
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where C u = Gsq/^t = R^sl/h 2 ~ 10 3 -10 4 . At later times, the alternating 
magnetic field rapidly decays because of diffusion and reconnection. The 
resulting maximum magnetic field strength achieved at to is given by 

£max~5o|CJ 1/2 , (34) 

where B is the external magnetic field; the magnetic field reverses at a small 
radial scale As ~ SolC^I -1 / 2 ~ 0.1 kpc. The magnetic pitch angle at to is of 
the order of \p\ ~ jC^I -1 / 2 < 1°, i.e., much smaller than the observed one. 
This picture cannot be reconciled with observations (cf. Kulsrud, 1999). It 
can be argued that streaming motions could make magnetic lines more open 
and parallel to the galactic spiral arms. However, then magnetic field will 
reverse on a small scale not only along radius, but also along azimuth. Such 
magnetic structures are quite different from what is observed. The moderate 
magnetic pitch angles observed in spiral galaxies are a direct indication that 
the regular magnetic field is not frozen into the interstellar gas and has to 
be maintained by the dynamo (Beck, 2000). 

5.2 The even (quadrupole) symmetry of magnetic field 
in the Milky Way 

One of the most convincing arguments in favour of the galactic dynamo the- 
ory comes from the symmetry of the observed regular magnetic field with 
respect to the Galactic equator in the Milky Way. The direction of the 
magnetic field is determined from Faraday rotation measures of the cosmic 
sources of polarized emission, pulsars and extragalactic radio sources. Since 
the Galactic magnetic field has a significant random component and extra- 
galactic radio sources can have their own (intrinsic) Faraday rotation, any 
meaningful conclusions about the Galactic magnetic field must rely on sta- 
tistically significant samples of Faraday rotation measures. Even though 
the quadrupole symmetry of the galactic magnetic fields has been widely 
accepted as a firmly established fact since mid-1970's, its objective obser- 
vational verification has been obtained only recently. The main problem 
here is that it is difficult to separate local (small-scale) and global magnetic 
structures in the observed picture. However, wavelet analysis of the Fara- 
day rotation measures of extragalactic radio sources has definitely confirmed 
that the horizontal components of the local regular magnetic field have even 
parity being similarly directed on both sides of the midplane (Frick et ai, 
2001, see Fig. lo~T| . 
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Figure 6: Upper panel: Faraday rotation measures of 551 extragalactic radio 
sources from the catalogue of Simard-Normandin & Kronberg (1980) shown 
in the (I, 6)-plane, where (I, b) are the Galactic longitude and latitude in a 
reference frame centred at the Sun with the Galactic center in the direction 
I = and Galactic midplane at b = 0. Positive (negative) RMs are shown 
with red (blue) circles whose radius indicates |RM| (radm -2 ) as shown to the 
right of the panel. The lower two panels show the wavelet transform of these 
data at scales 76° (middle panel) and 35° (lower panel) (Frick et ai, 2001). 
The transform at 76° has been obtained with the region of the Radio Loop I 
removed (this radio feature is a nearby supernova remnant). The wavelet 
transform at the scale 35° is dominated by local magneto-ionic features. 



The quadrupole symmetry is naturally explained by dynamo theory where 
even parity is strongly favoured against odd parity because the even field has 
twice larger scale in the vertical coordinate (see Sect. 14. Oj) . 

Primordial magnetic field twisted by differential rotation can have even 
vertical symmetry if it is parallel to the disc plane. However, then the field is 
rapidly destroyed by twisting and reconnection as described in Sect. 15.11 If, 
otherwise, the primordial field is parallel to the rotation axis and amplified by 
the vertical rotational shear dQ/dz (which, however, is insignificant within 
galactic discs, \z\ < h), it can avoid catastrophic decay (§3.11 in Moffatt, 
1978), but then it will have odd parity in z, which is ruled out by the observed 
parity of the Milky Way field. 

The derivation of the regular magnetic field of the Milky Way from Fara- 
day rotation measures of pulsars and extragalactic radio sources, RM, is 
complicated by the contribution of local magnetic perturbations, so it is dif- 
ficult to decide which features of the RM sky are due to the regular magnetic 
field and which are produced by localized magneto-ionic perturbations (e.g., 
supernova remnants). Therefore, the same observational data have lead dif- 
ferent authors to different conclusions (see Frick et al, 2001, for a recent 
review). Odd parity of the Galactic magnetic field has been suggested by 
Andreassian (1980, 1982) and, for the inner Galaxy, by Han et al. (1997). 
Quantitative methods of analysis (as opposed to the 'naked-eye' fitting of 
more or less arbitrarily selected models) are especially appropriate in this 
case. 

Unfortunately, it is difficult to determine the parity of magnetic field in 
external galaxies. In galaxies seen edge-on, the disc is depolarized, whereas 
Faraday rotation in the halo is weak. Beck et al. (1994) found weak evidence 
of even magnetic parity in the lower halo of NGC 253. The arrangement of 
polarization planes in the halo of NGC 4631 (Beck, 2000) is very suggestive 
of odd parity, but this does not exclude even parity in the disc. In galax- 
ies inclined to the line of sight, the amount of Faraday rotation produced 
by an odd (antisymmetric) magnetic field differs from zero because Faraday 
rotation and emission occur in the same volume; as a result, emission origi- 
nating at the far half of the galactic layer will have small or zero net rotation, 
whereas emission from the near half will have significant rotation produced 
by the unidirectional magnetic field in that half. Therefore, Faraday rota- 
tion measures produced by even and odd magnetic structures of the same 
strength only differ by a factor of two (Krause et al, 1989a; Sokoloff et al, 
1998) and it is difficult to distinguish between the two possibilities. 
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An interesting method to determine the parity of magnetic field in an 
external galaxy has been suggested by Han et al. (1998). These authors note 
that the contribution of the galaxy to the RM of a background radio source 
will be equal to the intrinsic RM of the galaxy if the magnetic field has even 
parity For odd parity, the galaxy will not contribute to the RM of a back- 
ground source, whereas any intrinsic RM will remain. The implementation 
of the method requires either a statistically significant sample of background 
sources or a single extended background source. 

5.3 The azimuthal structure 

Non-axisymmetric magnetic fields in a differentially rotating object are sub- 
ject to twisting and enhanced dissipation as described in Sect. 15.11 The 
dynamo can compensate for the losses, but axisymmetric magnetic fields are 
still easier to maintain (Radler, 1986). A few lowest non-axisymmetric modes 
with azimuthal wave numbers 



can be maintained in thin galactic disks where h <C So (§VTI.8 in Ruzmaikin 
et al, 1988). The WKB solution of the galactic acu-dynamo equations by 
Starchenko & Shukurov (1989) shows that the bisymmetric mode (m = 1) 
can grow provided 



which seems to be the case in some galaxies. These results indicate that it 
is natural to expect significant deviations from axial symmetry in magnetic 
fields of many spiral galaxies. However, the dominance of non-axisymmetric 
modes in most galaxies would be difficult to explain because the axisymmetric 
mode has the largest growth rate under typical conditions. 

Early interpretations of Faraday rotation in spiral galaxies were in strik- 
ing contrast with this picture, indicating strong dominance of bisymmetric 
magnetic structures (m — 1), B oc exp(z</>) with </> the azimuthal angle (Sofue 
et al., 1986), and this was considered to be a severe difficulty of the dynamo 
theory and an evidence of the primordial origin of galactic magnetic fields. It 
was suggested by Ruzmaikin et al. (1986) (see also Sawa & Fujimoto, 1986; 
Baryshnikova et al, 1987) that the bisymmetric magnetic structures can be 
interpreted as the m = 1 dynamo mode. However, despite effort, dynamo 
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models could not explain the apparent widespread dominance of bisymmetric 
magnetic structures. Paradoxically, what seemed to be a difficulty of the dy- 
namo theory has turned out to be its advantage as observations with better 
sensitivity and resolution and better interpretations have led to a dramatic 
revision of the observational picture. The present-day understanding is that 
modestly distorted axisymmetric magnetic structures occur in most galax- 
ies, wherein the dominant axisymmetric mode is mixed with weaker higher 
azimuthal modes (Beck et al, 1996; Beck, 2000). Among nearby galaxies, 
only M81 remains a candidate for a dominant bisymmetric magnetic struc- 
ture, but the data are old and this result needs to be reconsidered (Krause et 
al, 1989b); the interesting case of M51 is discussed below. Deviations from 
precise axial symmetry can result from the spiral pattern, asymmetry of the 
parent galaxy, etc. Dominant bisymmetric magnetic fields can be maintained 
by the dynamo action near the corotation radius due to a linear resonance 
with the spiral pattern (Mestel & Subramanian, 1991; Subramanian & Mes- 
tel, 1993; Moss, 1996) or nonlinear trapping of the field by the spiral pattern 
(Bykov et al, 1997). 

Twisting of a horizontal magnetic field by differential rotation generally 
produces a bisymmetric magnetic field, m — 1. Twisting of a horizontal 
primordial magnetic field can also produce an axisymmetric configuration 
near the galactic centre if the initial state is asymmetric (Sofue et al, 1986; 
Nordlund & Rognvaldsson, 2002), with a maximum of the primordial field 
displaced from the disc's rotation axis where the gas density is normally max- 
imum. Thus, the maximum of the primordial field required by this scenario 
has to occur at a different position than the maximum in the gas density. 
This can only occur if the primordial field is not frozen into the gas — other- 
wise the field strength scales as a positive power of gas density. The fact that 
magnetic fields in most spiral galaxies are nearly axisymmetric within large 
radius (in fact, in the whole galaxy) would require that this strong asymme- 
try in the initial state occurs systematically for all the galaxies, which would 
be difficult to explain. 

5.4 A composite magnetic structure in M51 and mag- 
netic reversals in the Milky Way 

A striking example of a complicated magnetic structure that can hardly be 
explained by any mechanism other than the dynamo has been revealed in 
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Figure 7: Left panels: The 

global magnetic structure of the 
galaxy M51, in the disc (upper 
left) and halo (bottom). Ar- 
rows show the direction and 
strength of the regular magnetic 
field on a polar grid shown su- 
perimposed on the optical im- 
age (Berkhuijsen et ai, 1997). 
The grid radii are 3, 6, 9, 
12 and 15kpc. Upper right 
panel: Magnetic field strength 
from the dynamo model for the 
disc of M51 (Bykov et al, 1997) 
is shown with shades of grey 
(darker shade means stronger 
field) . Magnetic field is reversed 
within the zero-level contour 
shown dashed; scale is given in 
kpc. The magnetic structure ro- 
tates rigidly together with the 
spiral pattern visible in the 
shades of grey. 
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the galaxy M51 by Berkhuijsen et al. (1997). These authors used radio 
polarization observations of the galaxy at wavelengths 2.8, 6.2, 18.0 and 
20.5cm (smoothed to a resolution of 3.5 kpc). The disc of this galaxy is 
not transparent to polarized radio emission at the two longer wavelengths. 
Therefore, it was possible to determine the magnetic field structure separately 
in two regions along the line the sight, which can be identified with the disc 
and halo of the galaxy. As shown in Fig. the regular magnetic fields in the 
disc is reversed in a region about 3 by 8 kpc in size extended along azimuth 
at galactocentric radii s = 3-6 kpc and azimuthal angles 300°-0 (shown with 
red arrows). A significant deviation from axial symmetry in the disc has been 
detected out to s = 9 kpc (in the azimuth range 160-260°), although it is too 
weak to result in a magnetic field reversal. The field reversal occurs around 
the corotation radius in M51, s ~ 6 kpc (i.e., the radius where the angular 
velocity of the spiral pattern is equal to that of the gas). 

A nonlinear dynamo model for M51 was developed by Bykov et al. (1997) 
who used the rotation curve of M51, with the pitch angle of the spiral arms 
— 15° and corotation radius 6 kpc. Figure [7| shows one of their solutions where 
a region with reversed magnetic field persists in the disc near the corotation 
radius of the spiral pattern. Near the corotation, a non-axisymmetric (bisym- 
metric) magnetic field can be trapped by the spiral pattern and maintained 
over the galactic lifetime. The effect is favoured by a smaller pitch angle of 
the spiral arms, thinner gaseous disc, weaker rotational shear and stronger 
spiral pattern. This nonaxisymmetric structure is arguably similar to the 
structure observed in M51. 

The regular magnetic field in the halo of M51 has a structure very different 
from that in the disc — the halo field is nearly axisymmetric and even directed 
oppositely to that in the disc in most of the galaxy. An external magnetic field 
should have a rather peculiar form to be twisted into such a configuration! 

Distinct azimuthal magnetic structures in the disc and the halo can be 
readily explained by dynamo theory as non-axisymmetric magnetic fields can 
be maintained only in the thin disc but not in the quasi-spherical halo where 
h ~ s and l-R^I 3> 1 in Eq. (|33|). Moreover, dynamo action in the disc and 
the halo can proceed almost independently of each other producing distinctly 
directed magnetic fields (Sokoloff & Shukurov, 1990). 

Another case of a regular magnetic field with unusual structure is the 
Milky Way where magnetic field reversals are observed along the galactocen- 
tric radius in the inner Galaxy between the Orion and Sagittarius arms at 
s ~ 7.9 kpc and, possibly, in the outer Galaxy between the Orion and Perseus 
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arms at s ~ 10.5 kpc (§3.8.2 in Beck et al, 1996, and Frick et al., 2001); see, 
however, Brown & Taylor 2001). The reversals were first interpreted as an 
indication of a global bisymmetric magnetic structure (Sofue & Fujimoto, 
1983), but it has been shown that dynamo-generated axisymmetric magnetic 
field can have reversals at the appropriate scale (Ruzmaikin et al., 1985; 
Poezd et al, 1993). Both interpretations presume that the reversals are of a 
global nature, i.e., they extend over the whole Galaxy to all azimuthal angles 
(or radii in the case of the bisymmetric structure). This leads to a question 
why reversals at this radial scale are not observed in any other galaxy (Beck, 
2000). Poezd et al. (1993) argue that the lifetime of the reversals is sensitive 
to subtle features of the rotation curve and the geometry of the ionized gas 
layer (see also Belyanin et al, 1994) and demonstrate that they are more 
probable to survive in the Milky Way than in, e.g., M31. 

However, the observational evidence of the reversals is restricted to a 
relatively small neighbourhood of the Sun, of at most 3-5 kpc along azimuth. 
It is therefore quite possible that the reversals are local and arise from a 
magnetic structure similar to that in the disc of M51 as shown in Fig. [7| 
The reversed field in the Solar neighbourhood has the same radial extent of 
2-3 kpc as in M51 and also occurs near the corotation radius. This possibility 
has not yet been explored; its observational verification would require careful 
analysis of pulsar Faraday rotation measures. 

5.5 The radial magnetic structure in M31 

An important clue to the origin of galactic magnetic field is provided by the 
magnetic ring in M31 (Beck, 1982), predicted by dynamo theory (Ruzmaikin 
& Shukurov, 1981). Both the large-scale magnetic field and the gas density 
in this galaxy have a maximum in the same annulus 8 < s < 12 kpc, with the 
apparent enhancement in the magnetic field strength by about 30% (Fletcher 
et al, 2004). The kinematic dynamo model of Ruzmaikin & Shukurov (1981) 
was based on the double-peaked rotation curve of shown in Fig. where 
rotational shear is strongly reduced at s = 2-6 kpc. As a result, R u is small 
and even positive in this radial range, so \D\ < \D CT \ and the dynamo cannot 
maintain any regular magnetic field at s — 2-6 kpc. 

An attractive aspect of this theory is that both magnetic and gas rings 
are attributed to the same feature of the rotation curve. Angular mo- 
mentum transport by viscous stress leads to matter inflow at a rate M = 
27r£z/T<9hi^V<9hi'S ~ 0.1M Q yr _1 , where ~ ?7t is the turbulent viscosity, 
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resulting in the radial inflow at a speed u s = M/27rsS with £ the gas surface 
density. In the nearly-rigidly rotating parts, u s is reduced and matter piles 
up outside this region producing gas ring. Gravitational torques from spiral 
arms can further enhance the inflow (see Moss et ai, 2000, for a discussion), 
so the total radial velocity is expected to be u s ~ lkms -1 at s = lOkpc. 

The double-peaked rotation curve of M31 is consistent with the existence 
of both magnetic and gas rings. The situation is different with the more 
recent rotation curve of Braun (1991) which does not have a double-peaked 
shape (Fig. EJ). The difference between the two rotation curves arises mainly 
from the fact that Braun allows for significant displacements of spiral arm 
segments from the galactic midplane: this results in a revision of the seg- 
ments' galactocentric distances for regions away from the major axis. We 
note that the CO velocity field at the major axis (Loinard et al, 1995) is 
compatible with a double-peaked rotation curve. 

With Braun's rotation curve, the magnetic field can concentrate into a 
ring mainly because the gas is in the ring and B oc p 1 / 2 as shown in Eq. ([24)1 . 
The dynamo model of Moss et al. (1998) based on the rotation curve of Braun 
(1991) has difficulties in reproducing a magnetic ring as well pronounced as 
implied by the observed amount of Faraday rotation. This has lead to an 
idea that magnetic field can be significant at s — 2-6 kpc in M31. This 
has prompted Han et al. (1998) to search for magnetic fields at s = 2-6 kpc 
that could have escaped detection because of reduced density of cosmic ray 
electrons at those radii. These authors have found that two out of three 
background polarized radio sources seen through that region of M31 have 
Faraday rotation measures compatible with the results of Moss et al. (1998). 
They further conclude that this indicates an even symmetry of the regular 
magnetic field. This is encouraging, but a statistically representative sample 
of background sources has to be used to reach definite conclusions because 
of their unknown intrinsic RM. 

With a double-peaked rotation curve, a primordial magnetic field with 
a uniform radial component could have been twisted to produce a magnetic 
ring by virtue of Eq. (|34j) . In this case the primordial and dynamo theories 
have similar problems and possibilities regarding the magnetic ring in M31. 

Lou & Fan (2000) attribute the magnetic ring in M31 to an axisymmetric 
mode of MHD density waves. Because of the axial symmetry of the wave, the 
magnetic field in the ring must be purely azimuthal, p = 0, in contrast to the 
observed structure with a significant pitch angle (Fig. EJ). Furthermore, the 
ring can hardly represent a wave packet as envisaged in this theory because 



44 



then its group velocity must be comparable to the Alfven speed of SOkms" 1 
(Lou & Fan, 1998) and so the ring should be travelling at this speed along 
radius to traverse 30kpc in 10 9 yr, a distance much larger than the ring 
radius. The implication would be that the ring is a transient with a short 
lifetime of order 3 x 10 8 yr. And, of course, the theory cannot explain the 
origin of an azimuthal magnetic field required to launch the wave packet. 

5.6 Strength of the regular magnetic field 

Interstellar regular magnetic fields are observed to be close to the energy 
equipartition with interstellar turbulence. This directly indicates that the 
regular magnetic field is coupled to the turbulent gas motions. (Note that 
£fl does not differ much from the turbulent velocity uq in Eq. 1)240.) To 
appreciate the importance of this conclusion, consider primordial magnetic 
field twisted by differential rotation. Its maximum strength given by Eq. (}3"4*|) 
as -B max ~ 10 2 -Bo is controlled by the strength of the primordial field B Q , 
and so this theory, if applicable, would result in stringent constraints on 
extragalactic magnetic fields. 

A striking evidence of the nontrivial behaviour of the large-scale galactic 
magnetic field is the so-called magnetic arms, first discovered by Beck & 
Hoernes (1996) in the nearby galaxy NGC 6946. We show in Fig. |H]a map of 
polarized radio emission from this galaxy (a tracer of the large-scale magnetic 
field strength) superimposed on the map in the Ha spectral line (a tracer of 
ionized gas). It is evident that magnetic field is stronger between the spiral 
arms of this galaxy, i.e., where the gas density (both total and ionized) is 
lower. This behaviour is just opposite to what is expected of a frozen-in 
magnetic field that scales with a power of gas density. The phenomenon of 
magnetic arms confirms in a spectacular manner that the large-scale magnetic 
field is not frozen into the interstellar gas, and therefore cannot be primordial. 
Shukurov (1998), Moss (1998) and Rohde et al. (1999) suggest an explanation 
of the magnetic arms in terms of the mean-field dynamo theory. In brief, they 
argue that dynamo number can be larger between the gaseous spiral arms, 
resulting in stronger dynamo action. 

The theory of MHD density waves relates magnetic field excess AB in 
spiral arms to the enhancement in stellar density, AB/ (B) = AS/ (E) (Lou & 
Fan, 1998), where S is the stellar surface density, AS is its excess in the spiral 
arms, and angular brackets denote azimuthal averaging. Arm intensities in 
magnetic field and stellar surface density in NGC 6946 have been estimated 
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Figure 8: Left panel: Magnetic arms in the galaxy NGC 6946: polarized 
intensity at the wavelength A = 6 cm (blue contours), a tracer of the large- 
scale magnetic field B, superimposed on the galactic image in the Ha spectral 
line of ionized hydrogen (grey scale). Red dashes indicate the orientation of 
the E'-vector of the polarized emission (parallel to the direction of intrinsic 
magnetic field if Faraday rotation is negligible), with length proportional to 
the fractional polarization — see Eq. (JSJ. The spiral arms visualized by Ha 
are the sites where gas density is maximum. The large-scale magnetic field 
is evidently stronger between the arms where gas density is lower. Right 
panel: As in the left panel, but now for the total synchrotron intensity, a 
tracer of the total magnetic field B 2 = B +b 2 . The total field is enhanced in 
the gaseous arms. Given that the large-scale field concentrates between the 
arms, this means that the random field is significantly stronger in the arms, 
a distribution very different from that of the large-scale field. The size of the 
beam in the radio maps is shown in the bottom left of each frame. (Images 
courtesy of R. Beck, MPIfR, Bonn.) 
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by Frick et al. (2000) who applied wavelet transform techniques to radio 
polarization maps and to the galaxy image in broadband red light, a tracer 
of stellar mass density. Their results indicate that the mean relative intensity 
of magnetic spiral arms remains rather constant with galactocentric radius 
at a level of 0.3-0.6. On the contrary, the relative strength of the stellar 
arms systematically grows with radius from very small values in the inner 
galaxy to 0.3-0.7 at s = 5-6 kpc, and then decreases to remain at a level of 
0.1-0.3 out to s = 12 kpc. The distinct magnitudes and radial trends in the 
strengths of magnetic and stellar arms in NGC 6946 do not seem to support 
the idea that the magnetic arms are due to MHD density waves. 

6 Elliptical galaxies 

Elliptical galaxies do not rotate fast enough, so they are ellipsoidal sys- 
tems without prominent disc components. The stellar population of ellip- 
tical galaxies is old and the interstellar gas is dilute (Mathews & Brighenti, 
2003). Therefore, both relativistic and thermal electrons have low density, 
and any synchrotron emission and Faraday rotation can only be weak. Nev- 
ertheless, there are several lines of evidence, albeit mostly indirect, suggest- 
ing significant magnetic fields in ellipticals (Moss & Shukurov, 1996; Math- 
ews & Brighenti, 1997). The magnetic field should be random, producing 
unpolarized synchrotron emission and fluctuating Faraday rotation. The 
r.m.s. Faraday rotation measure attributable to the ISM of the ellipticals is 
(RM) = 5-100 radm- 2 . 

6.1 Turbulent interstellar gas in elliptical galaxies 

Interstellar gas in elliptical galaxies is observed via its X-ray emission. Type I 
supernovae (SNe) (and also stellar winds and random motions of stars) heat 
the gas to the observed temperatures T ~ 10 7 K. It is natural to expect 
that a fraction 5 of the energy is converted into turbulent motions of the 
gas. The turbulent scale £ rs 400 pc is given by the diameter of a SN as 
it reaches pressure balance with the ambient medium whose typical density 
is n ~ 10~ 3 cm~ 3 . The balance between energy injection and dissipation 
rates yields a turbulent velocity of uq ~ 20kms _1 for 5 = 0.1, assuming the 
energy dissipation time r ~ £/u as for the Kolmogorov turbulence. This 
estimate of u is compatible with the constraint u <50kms _1 resulting from 
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the observed X-ray luminosity. Another driver of turbulence is the random 
motions of stars. These generate random vortical motions at a smaller scale 
and velocity, ~ 3pc and ~ 3kms _1 , respectively (Moss & Shukurov, 
1996). 

The driving force produced by an expanding quasi-spherical SN remnant 
is potential. The above estimates assume that the motions driven by the SNe 
are vortical, so r = £/u Q applies. In spiral galaxies, the potential (acoustic) 
motions are efficiently converted into vortical turbulence mainly due to the 
inhomogeneity of the ISM. The ISM in elliptical galaxies is hot and, presum- 
ably, rather homogeneous at kpc scales. Therefore, Moss & Shukurov (1996) 
suggested that SNe will drive sound-wave turbulence whose correlation time 
t is (u /c s )~ 2 £/c s ~ 3 x 10 7 yr rather than £/u , where c s ~ 300kms _1 is 
the speed of sound. However, Mathews & Brighenti (1997) noted that sound 
waves quickly dissipate, and so cannot form a pervasive turbulent velocity 
field. The nature of turbulence in elliptical galaxies needs to be studied 
further. 

6.2 The fluctuation dynamo in elliptical galaxies 

As in most astrophysical objects, Rm in elliptical galaxies by far exceeds 100, 
so fluctuation dynamo action in ellipticals is quite plausible (see Sect. 14.2(1 . 
The e-folding time of the random field in a vortical random flow is of the 
order of the eddy turnover time, r = 2 x 10 7 yr. The magnetic field is 
concentrated into flux ropes whose length and thickness are of order £ ~ 
400 pc and £Rm~ 1//2 ~ 40 pc. In the ropes, magnetic field is plausibly in 
equipartition with the turbulent kinetic energy, b ~ 0.3 fiG. 

Moss & Shukurov (1996) discuss a two-stage dynamo action by smaller 
scale vortical turbulence driven by random motions of the stars and, at larger 
scales, by the acoustic turbulence. However, the very existence of the acoustic 
turbulence in elliptical galaxies is questionable (Mathews & Brighenti, 1997). 

Faraday rotation measure produced within a single turbulent cell with the 
above parameters is RM — 0.81bn c £, so the net Faraday rotation from an en- 
semble of turbulent cells, observed at a resolution D such that D 3> £, is given 
by (RM) ~ RM oy jN/N D ~ RUq^L£/D , where N = Ljl is the number of 
cells along the path length L and No = (D/£) 2 is the number of cells in the 
resolution element. Thus, RMo — O.lradm -2 and (RM) ~ lradm -2 . Fara- 
day rotation can be stronger in the central regions where (RM) ~ 5 rad m -2 
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at a distance 8kpc from the galactic centre. These estimates agree fairly 
with the available observations. 

Magnetic field generation in elliptical galaxies was discussed by Lesch & 
Bender (1990), but they considered a mean-field dynamo that needs overall 
rotation which is not present in elliptical galaxies. The fluctuation dynamo in 
elliptical (radio) galaxies was simulated by De Young (1980), but these sim- 
ulations apparently had Rm < Rm cr as they resemble transient amplification 
of magnetic field by velocity shear rather than genuine dynamo action. 

7 Conclusions 

The observational picture of galactic magnetic fields is compatible with the 
mean-field dynamo theory in its simplest, quasi-kinematic form. It is im- 
portant to note that there is not much freedom in varying parameters of 
galactic dynamos as observations constrain them fairly tightly. Therefore, 
this agreement is not a result of a free-hand parameter adjustments. More- 
over, galactic dynamo theory has demonstrated its predictive power. For 
example, it has been clear to dynamo theorists since the early 1970's that 
the partially ionized Galactic disc has the scale height 0.4-0.5 kpc, i.e., signif- 
icantly larger than that of the neutral hydrogen layer (see §VI.2 in Ruzmaikin 
et ai, 1988, for a review), but the existence of this component of the inter- 
stellar medium was accepted by a broader astrophysical community only 15 
years later (Lockman, 1984). 

The agreement of the mean-field dynamo theory with observations dis- 
cussed in Sect. E] cannot be considered as a proof of its correctness — history 
of physics is familiar with concepts, such as the ether, that have proved to 
be irrelevant despite their perfect agreement with numerous experimental 
facts, before a single experiment refuted them. Nevertheless, the spectacular 
success of the dynamo theory when applied to galaxies warrants its careful 
treatment when compared with other theories. In particular, any rival theory 
has to be able to explain at least the same set of observational data as the 
dynamo theory. In this sense, there are no fair rivals to the dynamo theory. 

The main current difficulty of the galactic dynamo theory is our lack of 
understanding of its nonlinear form. It is important to avoid an unjustified 
extension to real galaxies of results obtained for highly idealized systems. In 
particular, the disc of a spiral galaxy is not a closed system. The significance 
of the disc-halo connection for the mean-field galactic dynamos, touched upon 
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in Sect. 14 .3\ should be carefully investigated. 

Another outstanding problem in theory and observations of galactic mag- 
netic fields is the effect of the multi-phase structure of the interstellar medium 
on the magnetic field. The effect of magnetic fields on the multi-phase struc- 
ture is also expected to be very significant, and also poorly understood. The 
strength of the magnetic field in the hot phase is not known. The detailed na- 
ture of the balance between cosmic rays and magnetic field has to be clarified: 
it is unclear whether this balance is maintained pointwise (at each location) 
or only on average (e.g., at scales exceeding the diffusion scale of cosmic ray 
particles). The answer to this question is essential for the interpretation of 
the synchrotron emission from spiral galaxies. Progress in this direction will 
eventually make theory of galactic magnetic fields an integral part of galactic 
astrophysics. 
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